3.5.18 \(\int \frac {\cos ^3(c+d x) (A+B \cos (c+d x)+C \cos ^2(c+d x))}{(a+a \cos (c+d x))^{5/2}} \, dx\) [418]

3.5.18.1 Optimal result
3.5.18.2 Mathematica [A] (verified)
3.5.18.3 Rubi [A] (verified)
3.5.18.4 Maple [B] (verified)
3.5.18.5 Fricas [A] (verification not implemented)
3.5.18.6 Sympy [F(-1)]
3.5.18.7 Maxima [F(-1)]
3.5.18.8 Giac [A] (verification not implemented)
3.5.18.9 Mupad [F(-1)]

3.5.18.1 Optimal result

Integrand size = 43, antiderivative size = 277 \[ \int \frac {\cos ^3(c+d x) \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right )}{(a+a \cos (c+d x))^{5/2}} \, dx=-\frac {(75 A-163 B+283 C) \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {a+a \cos (c+d x)}}\right )}{16 \sqrt {2} a^{5/2} d}-\frac {(A-B+C) \cos ^4(c+d x) \sin (c+d x)}{4 d (a+a \cos (c+d x))^{5/2}}-\frac {(5 A-13 B+21 C) \cos ^3(c+d x) \sin (c+d x)}{16 a d (a+a \cos (c+d x))^{3/2}}+\frac {(465 A-985 B+1729 C) \sin (c+d x)}{120 a^2 d \sqrt {a+a \cos (c+d x)}}+\frac {(45 A-85 B+157 C) \cos ^2(c+d x) \sin (c+d x)}{80 a^2 d \sqrt {a+a \cos (c+d x)}}-\frac {(195 A-475 B+787 C) \sqrt {a+a \cos (c+d x)} \sin (c+d x)}{240 a^3 d} \]

output
-1/4*(A-B+C)*cos(d*x+c)^4*sin(d*x+c)/d/(a+a*cos(d*x+c))^(5/2)-1/16*(5*A-13 
*B+21*C)*cos(d*x+c)^3*sin(d*x+c)/a/d/(a+a*cos(d*x+c))^(3/2)-1/32*(75*A-163 
*B+283*C)*arctanh(1/2*sin(d*x+c)*a^(1/2)*2^(1/2)/(a+a*cos(d*x+c))^(1/2))/a 
^(5/2)/d*2^(1/2)+1/120*(465*A-985*B+1729*C)*sin(d*x+c)/a^2/d/(a+a*cos(d*x+ 
c))^(1/2)+1/80*(45*A-85*B+157*C)*cos(d*x+c)^2*sin(d*x+c)/a^2/d/(a+a*cos(d* 
x+c))^(1/2)-1/240*(195*A-475*B+787*C)*sin(d*x+c)*(a+a*cos(d*x+c))^(1/2)/a^ 
3/d
 
3.5.18.2 Mathematica [A] (verified)

Time = 1.70 (sec) , antiderivative size = 152, normalized size of antiderivative = 0.55 \[ \int \frac {\cos ^3(c+d x) \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right )}{(a+a \cos (c+d x))^{5/2}} \, dx=\frac {-30 (75 A-163 B+283 C) \text {arctanh}\left (\sin \left (\frac {1}{2} (c+d x)\right )\right ) \cos ^3\left (\frac {1}{2} (c+d x)\right )+(975 A-1895 B+3491 C+5 (255 A-479 B+887 C) \cos (c+d x)+16 (15 A-25 B+52 C) \cos (2 (c+d x))+40 B \cos (3 (c+d x))-40 C \cos (3 (c+d x))+12 C \cos (4 (c+d x))) \tan \left (\frac {1}{2} (c+d x)\right )}{240 a d (a (1+\cos (c+d x)))^{3/2}} \]

input
Integrate[(Cos[c + d*x]^3*(A + B*Cos[c + d*x] + C*Cos[c + d*x]^2))/(a + a* 
Cos[c + d*x])^(5/2),x]
 
output
(-30*(75*A - 163*B + 283*C)*ArcTanh[Sin[(c + d*x)/2]]*Cos[(c + d*x)/2]^3 + 
 (975*A - 1895*B + 3491*C + 5*(255*A - 479*B + 887*C)*Cos[c + d*x] + 16*(1 
5*A - 25*B + 52*C)*Cos[2*(c + d*x)] + 40*B*Cos[3*(c + d*x)] - 40*C*Cos[3*( 
c + d*x)] + 12*C*Cos[4*(c + d*x)])*Tan[(c + d*x)/2])/(240*a*d*(a*(1 + Cos[ 
c + d*x]))^(3/2))
 
3.5.18.3 Rubi [A] (verified)

Time = 1.84 (sec) , antiderivative size = 303, normalized size of antiderivative = 1.09, number of steps used = 20, number of rules used = 19, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.442, Rules used = {3042, 3520, 27, 3042, 3456, 27, 3042, 3462, 27, 3042, 3447, 3042, 3502, 27, 3042, 3230, 3042, 3128, 219}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\cos ^3(c+d x) \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right )}{(a \cos (c+d x)+a)^{5/2}} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {\sin \left (c+d x+\frac {\pi }{2}\right )^3 \left (A+B \sin \left (c+d x+\frac {\pi }{2}\right )+C \sin \left (c+d x+\frac {\pi }{2}\right )^2\right )}{\left (a \sin \left (c+d x+\frac {\pi }{2}\right )+a\right )^{5/2}}dx\)

\(\Big \downarrow \) 3520

\(\displaystyle \frac {\int \frac {\cos ^3(c+d x) (8 a (B-C)+a (5 A-5 B+13 C) \cos (c+d x))}{2 (\cos (c+d x) a+a)^{3/2}}dx}{4 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\int \frac {\cos ^3(c+d x) (8 a (B-C)+a (5 A-5 B+13 C) \cos (c+d x))}{(\cos (c+d x) a+a)^{3/2}}dx}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\int \frac {\sin \left (c+d x+\frac {\pi }{2}\right )^3 \left (8 a (B-C)+a (5 A-5 B+13 C) \sin \left (c+d x+\frac {\pi }{2}\right )\right )}{\left (\sin \left (c+d x+\frac {\pi }{2}\right ) a+a\right )^{3/2}}dx}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3456

\(\displaystyle \frac {\frac {\int -\frac {\cos ^2(c+d x) \left (6 a^2 (5 A-13 B+21 C)-a^2 (45 A-85 B+157 C) \cos (c+d x)\right )}{2 \sqrt {\cos (c+d x) a+a}}dx}{2 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {-\frac {\int \frac {\cos ^2(c+d x) \left (6 a^2 (5 A-13 B+21 C)-a^2 (45 A-85 B+157 C) \cos (c+d x)\right )}{\sqrt {\cos (c+d x) a+a}}dx}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {-\frac {\int \frac {\sin \left (c+d x+\frac {\pi }{2}\right )^2 \left (6 a^2 (5 A-13 B+21 C)-a^2 (45 A-85 B+157 C) \sin \left (c+d x+\frac {\pi }{2}\right )\right )}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}dx}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3462

\(\displaystyle \frac {-\frac {\frac {2 \int -\frac {\cos (c+d x) \left (4 a^3 (45 A-85 B+157 C)-a^3 (195 A-475 B+787 C) \cos (c+d x)\right )}{2 \sqrt {\cos (c+d x) a+a}}dx}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {-\frac {-\frac {\int \frac {\cos (c+d x) \left (4 a^3 (45 A-85 B+157 C)-a^3 (195 A-475 B+787 C) \cos (c+d x)\right )}{\sqrt {\cos (c+d x) a+a}}dx}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {-\frac {-\frac {\int \frac {\sin \left (c+d x+\frac {\pi }{2}\right ) \left (4 a^3 (45 A-85 B+157 C)-a^3 (195 A-475 B+787 C) \sin \left (c+d x+\frac {\pi }{2}\right )\right )}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}dx}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3447

\(\displaystyle \frac {-\frac {-\frac {\int \frac {4 a^3 (45 A-85 B+157 C) \cos (c+d x)-a^3 (195 A-475 B+787 C) \cos ^2(c+d x)}{\sqrt {\cos (c+d x) a+a}}dx}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {-\frac {-\frac {\int \frac {4 a^3 (45 A-85 B+157 C) \sin \left (c+d x+\frac {\pi }{2}\right )-a^3 (195 A-475 B+787 C) \sin \left (c+d x+\frac {\pi }{2}\right )^2}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}dx}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3502

\(\displaystyle \frac {-\frac {-\frac {\frac {2 \int -\frac {a^4 (195 A-475 B+787 C)-2 a^4 (465 A-985 B+1729 C) \cos (c+d x)}{2 \sqrt {\cos (c+d x) a+a}}dx}{3 a}-\frac {2 a^2 (195 A-475 B+787 C) \sin (c+d x) \sqrt {a \cos (c+d x)+a}}{3 d}}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {-\frac {-\frac {-\frac {\int \frac {a^4 (195 A-475 B+787 C)-2 a^4 (465 A-985 B+1729 C) \cos (c+d x)}{\sqrt {\cos (c+d x) a+a}}dx}{3 a}-\frac {2 a^2 (195 A-475 B+787 C) \sin (c+d x) \sqrt {a \cos (c+d x)+a}}{3 d}}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {-\frac {-\frac {-\frac {\int \frac {a^4 (195 A-475 B+787 C)-2 a^4 (465 A-985 B+1729 C) \sin \left (c+d x+\frac {\pi }{2}\right )}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}dx}{3 a}-\frac {2 a^2 (195 A-475 B+787 C) \sin (c+d x) \sqrt {a \cos (c+d x)+a}}{3 d}}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3230

\(\displaystyle \frac {-\frac {-\frac {-\frac {15 a^4 (75 A-163 B+283 C) \int \frac {1}{\sqrt {\cos (c+d x) a+a}}dx-\frac {4 a^4 (465 A-985 B+1729 C) \sin (c+d x)}{d \sqrt {a \cos (c+d x)+a}}}{3 a}-\frac {2 a^2 (195 A-475 B+787 C) \sin (c+d x) \sqrt {a \cos (c+d x)+a}}{3 d}}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {-\frac {-\frac {-\frac {15 a^4 (75 A-163 B+283 C) \int \frac {1}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}dx-\frac {4 a^4 (465 A-985 B+1729 C) \sin (c+d x)}{d \sqrt {a \cos (c+d x)+a}}}{3 a}-\frac {2 a^2 (195 A-475 B+787 C) \sin (c+d x) \sqrt {a \cos (c+d x)+a}}{3 d}}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 3128

\(\displaystyle \frac {-\frac {-\frac {-\frac {-\frac {30 a^4 (75 A-163 B+283 C) \int \frac {1}{2 a-\frac {a^2 \sin ^2(c+d x)}{\cos (c+d x) a+a}}d\left (-\frac {a \sin (c+d x)}{\sqrt {\cos (c+d x) a+a}}\right )}{d}-\frac {4 a^4 (465 A-985 B+1729 C) \sin (c+d x)}{d \sqrt {a \cos (c+d x)+a}}}{3 a}-\frac {2 a^2 (195 A-475 B+787 C) \sin (c+d x) \sqrt {a \cos (c+d x)+a}}{3 d}}{5 a}-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {-\frac {-\frac {2 a^2 (45 A-85 B+157 C) \sin (c+d x) \cos ^2(c+d x)}{5 d \sqrt {a \cos (c+d x)+a}}-\frac {-\frac {2 a^2 (195 A-475 B+787 C) \sin (c+d x) \sqrt {a \cos (c+d x)+a}}{3 d}-\frac {\frac {15 \sqrt {2} a^{7/2} (75 A-163 B+283 C) \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {a \cos (c+d x)+a}}\right )}{d}-\frac {4 a^4 (465 A-985 B+1729 C) \sin (c+d x)}{d \sqrt {a \cos (c+d x)+a}}}{3 a}}{5 a}}{4 a^2}-\frac {a (5 A-13 B+21 C) \sin (c+d x) \cos ^3(c+d x)}{2 d (a \cos (c+d x)+a)^{3/2}}}{8 a^2}-\frac {(A-B+C) \sin (c+d x) \cos ^4(c+d x)}{4 d (a \cos (c+d x)+a)^{5/2}}\)

input
Int[(Cos[c + d*x]^3*(A + B*Cos[c + d*x] + C*Cos[c + d*x]^2))/(a + a*Cos[c 
+ d*x])^(5/2),x]
 
output
-1/4*((A - B + C)*Cos[c + d*x]^4*Sin[c + d*x])/(d*(a + a*Cos[c + d*x])^(5/ 
2)) + (-1/2*(a*(5*A - 13*B + 21*C)*Cos[c + d*x]^3*Sin[c + d*x])/(d*(a + a* 
Cos[c + d*x])^(3/2)) - ((-2*a^2*(45*A - 85*B + 157*C)*Cos[c + d*x]^2*Sin[c 
 + d*x])/(5*d*Sqrt[a + a*Cos[c + d*x]]) - ((-2*a^2*(195*A - 475*B + 787*C) 
*Sqrt[a + a*Cos[c + d*x]]*Sin[c + d*x])/(3*d) - ((15*Sqrt[2]*a^(7/2)*(75*A 
 - 163*B + 283*C)*ArcTanh[(Sqrt[a]*Sin[c + d*x])/(Sqrt[2]*Sqrt[a + a*Cos[c 
 + d*x]])])/d - (4*a^4*(465*A - 985*B + 1729*C)*Sin[c + d*x])/(d*Sqrt[a + 
a*Cos[c + d*x]]))/(3*a))/(5*a))/(4*a^2))/(8*a^2)
 

3.5.18.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 219
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))* 
ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && NegQ[a/b] && (Gt 
Q[a, 0] || LtQ[b, 0])
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3128
Int[1/Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[-2/d 
Subst[Int[1/(2*a - x^2), x], x, b*(Cos[c + d*x]/Sqrt[a + b*Sin[c + d*x]])], 
 x] /; FreeQ[{a, b, c, d}, x] && EqQ[a^2 - b^2, 0]
 

rule 3230
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + 
(f_.)*(x_)]), x_Symbol] :> Simp[(-d)*Cos[e + f*x]*((a + b*Sin[e + f*x])^m/( 
f*(m + 1))), x] + Simp[(a*d*m + b*c*(m + 1))/(b*(m + 1))   Int[(a + b*Sin[e 
 + f*x])^m, x], x] /; FreeQ[{a, b, c, d, e, f, m}, x] && NeQ[b*c - a*d, 0] 
&& EqQ[a^2 - b^2, 0] &&  !LtQ[m, -2^(-1)]
 

rule 3447
Int[((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((A_.) + (B_.)*sin[(e_.) 
+ (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Int[(a 
 + b*Sin[e + f*x])^m*(A*c + (B*c + A*d)*Sin[e + f*x] + B*d*Sin[e + f*x]^2), 
 x] /; FreeQ[{a, b, c, d, e, f, A, B, m}, x] && NeQ[b*c - a*d, 0]
 

rule 3456
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + 
(f_.)*(x_)])*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Sim 
p[(A*b - a*B)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^n/( 
a*f*(2*m + 1))), x] - Simp[1/(a*b*(2*m + 1))   Int[(a + b*Sin[e + f*x])^(m 
+ 1)*(c + d*Sin[e + f*x])^(n - 1)*Simp[A*(a*d*n - b*c*(m + 1)) - B*(a*c*m + 
 b*d*n) - d*(a*B*(m - n) + A*b*(m + n + 1))*Sin[e + f*x], x], x], x] /; Fre 
eQ[{a, b, c, d, e, f, A, B}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] & 
& NeQ[c^2 - d^2, 0] && LtQ[m, -2^(-1)] && GtQ[n, 0] && IntegerQ[2*m] && (In 
tegerQ[2*n] || EqQ[c, 0])
 

rule 3462
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + 
(f_.)*(x_)])*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Sim 
p[(-B)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^n/(f*(m + 
n + 1))), x] + Simp[1/(b*(m + n + 1))   Int[(a + b*Sin[e + f*x])^m*(c + d*S 
in[e + f*x])^(n - 1)*Simp[A*b*c*(m + n + 1) + B*(a*c*m + b*d*n) + (A*b*d*(m 
 + n + 1) + B*(a*d*m + b*c*n))*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, c, 
d, e, f, A, B, m}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 
- d^2, 0] && GtQ[n, 0] && (IntegerQ[n] || EqQ[m + 1/2, 0])
 

rule 3502
Int[((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((A_.) + (B_.)*sin[(e_.) 
+ (f_.)*(x_)] + (C_.)*sin[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> Simp[(-C)*Co 
s[e + f*x]*((a + b*Sin[e + f*x])^(m + 1)/(b*f*(m + 2))), x] + Simp[1/(b*(m 
+ 2))   Int[(a + b*Sin[e + f*x])^m*Simp[A*b*(m + 2) + b*C*(m + 1) + (b*B*(m 
 + 2) - a*C)*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, e, f, A, B, C, m}, x] 
 &&  !LtQ[m, -1]
 

rule 3520
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.) + 
(f_.)*(x_)])^(n_.)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)] + (C_.)*sin[(e_.) 
 + (f_.)*(x_)]^2), x_Symbol] :> Simp[(a*A - b*B + a*C)*Cos[e + f*x]*(a + b* 
Sin[e + f*x])^m*((c + d*Sin[e + f*x])^(n + 1)/(f*(b*c - a*d)*(2*m + 1))), x 
] + Simp[1/(b*(b*c - a*d)*(2*m + 1))   Int[(a + b*Sin[e + f*x])^(m + 1)*(c 
+ d*Sin[e + f*x])^n*Simp[A*(a*c*(m + 1) - b*d*(2*m + n + 2)) + B*(b*c*m + a 
*d*(n + 1)) - C*(a*c*m + b*d*(n + 1)) + (d*(a*A - b*B)*(m + n + 2) + C*(b*c 
*(2*m + 1) - a*d*(m - n - 1)))*Sin[e + f*x], x], x], x] /; FreeQ[{a, b, c, 
d, e, f, A, B, C, n}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c 
^2 - d^2, 0] && LtQ[m, -2^(-1)]
 
3.5.18.4 Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(616\) vs. \(2(246)=492\).

Time = 12.45 (sec) , antiderivative size = 617, normalized size of antiderivative = 2.23

method result size
default \(\frac {\sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, \left (768 C \sqrt {2}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, \sqrt {a}\, \left (\cos ^{8}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )+640 B \sqrt {2}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, \sqrt {a}\, \left (\cos ^{6}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )-2176 C \sqrt {a}\, \sqrt {2}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, \left (\cos ^{6}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )-1125 A \sqrt {2}\, \ln \left (\frac {4 \sqrt {a}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}+4 a}{\cos \left (\frac {d x}{2}+\frac {c}{2}\right )}\right ) \left (\cos ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right ) a +2445 B \sqrt {2}\, \ln \left (\frac {4 \sqrt {a}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}+4 a}{\cos \left (\frac {d x}{2}+\frac {c}{2}\right )}\right ) \left (\cos ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right ) a -4245 C \sqrt {2}\, \ln \left (\frac {4 \sqrt {a}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}+4 a}{\cos \left (\frac {d x}{2}+\frac {c}{2}\right )}\right ) \left (\cos ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right ) a +960 A \sqrt {2}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, \sqrt {a}\, \left (\cos ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )-2560 B \left (\cos ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right ) \sqrt {2}\, \sqrt {a}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}+5248 C \sqrt {a}\, \sqrt {2}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, \left (\cos ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )+315 A \left (\cos ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right ) \sqrt {2}\, \sqrt {a}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}-435 B \left (\cos ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right ) \sqrt {2}\, \sqrt {a}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}+555 C \sqrt {a}\, \sqrt {2}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, \left (\cos ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )-30 A \sqrt {a}\, \sqrt {2}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}+30 B \sqrt {2}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, \sqrt {a}-30 C \sqrt {2}\, \sqrt {a \left (\sin ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, \sqrt {a}\right )}{480 \cos \left (\frac {d x}{2}+\frac {c}{2}\right )^{3} a^{\frac {7}{2}} \sin \left (\frac {d x}{2}+\frac {c}{2}\right ) \sqrt {a \left (\cos ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}\, d}\) \(617\)
parts \(\text {Expression too large to display}\) \(728\)

input
int(cos(d*x+c)^3*(A+B*cos(d*x+c)+C*cos(d*x+c)^2)/(a+cos(d*x+c)*a)^(5/2),x, 
method=_RETURNVERBOSE)
 
output
1/480/cos(1/2*d*x+1/2*c)^3*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)*(768*C*2^(1/2)*( 
a*sin(1/2*d*x+1/2*c)^2)^(1/2)*a^(1/2)*cos(1/2*d*x+1/2*c)^8+640*B*2^(1/2)*( 
a*sin(1/2*d*x+1/2*c)^2)^(1/2)*a^(1/2)*cos(1/2*d*x+1/2*c)^6-2176*C*a^(1/2)* 
2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)*cos(1/2*d*x+1/2*c)^6-1125*A*2^(1/2) 
*ln(2*(2*a^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)+2*a)/cos(1/2*d*x+1/2*c))*c 
os(1/2*d*x+1/2*c)^4*a+2445*B*2^(1/2)*ln(2*(2*a^(1/2)*(a*sin(1/2*d*x+1/2*c) 
^2)^(1/2)+2*a)/cos(1/2*d*x+1/2*c))*cos(1/2*d*x+1/2*c)^4*a-4245*C*2^(1/2)*l 
n(2*(2*a^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)+2*a)/cos(1/2*d*x+1/2*c))*cos 
(1/2*d*x+1/2*c)^4*a+960*A*2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)*a^(1/2)*c 
os(1/2*d*x+1/2*c)^4-2560*B*cos(1/2*d*x+1/2*c)^4*2^(1/2)*a^(1/2)*(a*sin(1/2 
*d*x+1/2*c)^2)^(1/2)+5248*C*a^(1/2)*2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2) 
*cos(1/2*d*x+1/2*c)^4+315*A*cos(1/2*d*x+1/2*c)^2*2^(1/2)*a^(1/2)*(a*sin(1/ 
2*d*x+1/2*c)^2)^(1/2)-435*B*cos(1/2*d*x+1/2*c)^2*2^(1/2)*a^(1/2)*(a*sin(1/ 
2*d*x+1/2*c)^2)^(1/2)+555*C*a^(1/2)*2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2) 
*cos(1/2*d*x+1/2*c)^2-30*A*a^(1/2)*2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)+ 
30*B*2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)*a^(1/2)-30*C*2^(1/2)*(a*sin(1/ 
2*d*x+1/2*c)^2)^(1/2)*a^(1/2))/a^(7/2)/sin(1/2*d*x+1/2*c)/(a*cos(1/2*d*x+1 
/2*c)^2)^(1/2)/d
 
3.5.18.5 Fricas [A] (verification not implemented)

Time = 0.29 (sec) , antiderivative size = 291, normalized size of antiderivative = 1.05 \[ \int \frac {\cos ^3(c+d x) \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right )}{(a+a \cos (c+d x))^{5/2}} \, dx=\frac {15 \, \sqrt {2} {\left ({\left (75 \, A - 163 \, B + 283 \, C\right )} \cos \left (d x + c\right )^{3} + 3 \, {\left (75 \, A - 163 \, B + 283 \, C\right )} \cos \left (d x + c\right )^{2} + 3 \, {\left (75 \, A - 163 \, B + 283 \, C\right )} \cos \left (d x + c\right ) + 75 \, A - 163 \, B + 283 \, C\right )} \sqrt {a} \log \left (-\frac {a \cos \left (d x + c\right )^{2} + 2 \, \sqrt {2} \sqrt {a \cos \left (d x + c\right ) + a} \sqrt {a} \sin \left (d x + c\right ) - 2 \, a \cos \left (d x + c\right ) - 3 \, a}{\cos \left (d x + c\right )^{2} + 2 \, \cos \left (d x + c\right ) + 1}\right ) + 4 \, {\left (96 \, C \cos \left (d x + c\right )^{4} + 160 \, {\left (B - C\right )} \cos \left (d x + c\right )^{3} + 32 \, {\left (15 \, A - 25 \, B + 49 \, C\right )} \cos \left (d x + c\right )^{2} + 5 \, {\left (255 \, A - 503 \, B + 911 \, C\right )} \cos \left (d x + c\right ) + 735 \, A - 1495 \, B + 2671 \, C\right )} \sqrt {a \cos \left (d x + c\right ) + a} \sin \left (d x + c\right )}{960 \, {\left (a^{3} d \cos \left (d x + c\right )^{3} + 3 \, a^{3} d \cos \left (d x + c\right )^{2} + 3 \, a^{3} d \cos \left (d x + c\right ) + a^{3} d\right )}} \]

input
integrate(cos(d*x+c)^3*(A+B*cos(d*x+c)+C*cos(d*x+c)^2)/(a+a*cos(d*x+c))^(5 
/2),x, algorithm="fricas")
 
output
1/960*(15*sqrt(2)*((75*A - 163*B + 283*C)*cos(d*x + c)^3 + 3*(75*A - 163*B 
 + 283*C)*cos(d*x + c)^2 + 3*(75*A - 163*B + 283*C)*cos(d*x + c) + 75*A - 
163*B + 283*C)*sqrt(a)*log(-(a*cos(d*x + c)^2 + 2*sqrt(2)*sqrt(a*cos(d*x + 
 c) + a)*sqrt(a)*sin(d*x + c) - 2*a*cos(d*x + c) - 3*a)/(cos(d*x + c)^2 + 
2*cos(d*x + c) + 1)) + 4*(96*C*cos(d*x + c)^4 + 160*(B - C)*cos(d*x + c)^3 
 + 32*(15*A - 25*B + 49*C)*cos(d*x + c)^2 + 5*(255*A - 503*B + 911*C)*cos( 
d*x + c) + 735*A - 1495*B + 2671*C)*sqrt(a*cos(d*x + c) + a)*sin(d*x + c)) 
/(a^3*d*cos(d*x + c)^3 + 3*a^3*d*cos(d*x + c)^2 + 3*a^3*d*cos(d*x + c) + a 
^3*d)
 
3.5.18.6 Sympy [F(-1)]

Timed out. \[ \int \frac {\cos ^3(c+d x) \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right )}{(a+a \cos (c+d x))^{5/2}} \, dx=\text {Timed out} \]

input
integrate(cos(d*x+c)**3*(A+B*cos(d*x+c)+C*cos(d*x+c)**2)/(a+a*cos(d*x+c))* 
*(5/2),x)
 
output
Timed out
 
3.5.18.7 Maxima [F(-1)]

Timed out. \[ \int \frac {\cos ^3(c+d x) \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right )}{(a+a \cos (c+d x))^{5/2}} \, dx=\text {Timed out} \]

input
integrate(cos(d*x+c)^3*(A+B*cos(d*x+c)+C*cos(d*x+c)^2)/(a+a*cos(d*x+c))^(5 
/2),x, algorithm="maxima")
 
output
Timed out
 
3.5.18.8 Giac [A] (verification not implemented)

Time = 15.28 (sec) , antiderivative size = 374, normalized size of antiderivative = 1.35 \[ \int \frac {\cos ^3(c+d x) \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right )}{(a+a \cos (c+d x))^{5/2}} \, dx=-\frac {\frac {15 \, \sqrt {2} {\left (75 \, A \sqrt {a} - 163 \, B \sqrt {a} + 283 \, C \sqrt {a}\right )} \log \left (\sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) + 1\right )}{a^{3} \mathrm {sgn}\left (\cos \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )\right )} - \frac {15 \, \sqrt {2} {\left (75 \, A \sqrt {a} - 163 \, B \sqrt {a} + 283 \, C \sqrt {a}\right )} \log \left (-\sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) + 1\right )}{a^{3} \mathrm {sgn}\left (\cos \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )\right )} + \frac {30 \, {\left (21 \, \sqrt {2} A \sqrt {a} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{3} - 29 \, \sqrt {2} B \sqrt {a} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{3} + 37 \, \sqrt {2} C \sqrt {a} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{3} - 19 \, \sqrt {2} A \sqrt {a} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) + 27 \, \sqrt {2} B \sqrt {a} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) - 35 \, \sqrt {2} C \sqrt {a} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )\right )}}{{\left (\sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{2} - 1\right )}^{2} a^{3} \mathrm {sgn}\left (\cos \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )\right )} - \frac {128 \, \sqrt {2} {\left (12 \, C a^{\frac {25}{2}} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{5} - 10 \, B a^{\frac {25}{2}} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{3} + 10 \, C a^{\frac {25}{2}} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )^{3} + 15 \, A a^{\frac {25}{2}} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) - 30 \, B a^{\frac {25}{2}} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right ) + 60 \, C a^{\frac {25}{2}} \sin \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )\right )}}{a^{15} \mathrm {sgn}\left (\cos \left (\frac {1}{2} \, d x + \frac {1}{2} \, c\right )\right )}}{960 \, d} \]

input
integrate(cos(d*x+c)^3*(A+B*cos(d*x+c)+C*cos(d*x+c)^2)/(a+a*cos(d*x+c))^(5 
/2),x, algorithm="giac")
 
output
-1/960*(15*sqrt(2)*(75*A*sqrt(a) - 163*B*sqrt(a) + 283*C*sqrt(a))*log(sin( 
1/2*d*x + 1/2*c) + 1)/(a^3*sgn(cos(1/2*d*x + 1/2*c))) - 15*sqrt(2)*(75*A*s 
qrt(a) - 163*B*sqrt(a) + 283*C*sqrt(a))*log(-sin(1/2*d*x + 1/2*c) + 1)/(a^ 
3*sgn(cos(1/2*d*x + 1/2*c))) + 30*(21*sqrt(2)*A*sqrt(a)*sin(1/2*d*x + 1/2* 
c)^3 - 29*sqrt(2)*B*sqrt(a)*sin(1/2*d*x + 1/2*c)^3 + 37*sqrt(2)*C*sqrt(a)* 
sin(1/2*d*x + 1/2*c)^3 - 19*sqrt(2)*A*sqrt(a)*sin(1/2*d*x + 1/2*c) + 27*sq 
rt(2)*B*sqrt(a)*sin(1/2*d*x + 1/2*c) - 35*sqrt(2)*C*sqrt(a)*sin(1/2*d*x + 
1/2*c))/((sin(1/2*d*x + 1/2*c)^2 - 1)^2*a^3*sgn(cos(1/2*d*x + 1/2*c))) - 1 
28*sqrt(2)*(12*C*a^(25/2)*sin(1/2*d*x + 1/2*c)^5 - 10*B*a^(25/2)*sin(1/2*d 
*x + 1/2*c)^3 + 10*C*a^(25/2)*sin(1/2*d*x + 1/2*c)^3 + 15*A*a^(25/2)*sin(1 
/2*d*x + 1/2*c) - 30*B*a^(25/2)*sin(1/2*d*x + 1/2*c) + 60*C*a^(25/2)*sin(1 
/2*d*x + 1/2*c))/(a^15*sgn(cos(1/2*d*x + 1/2*c))))/d
 
3.5.18.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\cos ^3(c+d x) \left (A+B \cos (c+d x)+C \cos ^2(c+d x)\right )}{(a+a \cos (c+d x))^{5/2}} \, dx=\int \frac {{\cos \left (c+d\,x\right )}^3\,\left (C\,{\cos \left (c+d\,x\right )}^2+B\,\cos \left (c+d\,x\right )+A\right )}{{\left (a+a\,\cos \left (c+d\,x\right )\right )}^{5/2}} \,d x \]

input
int((cos(c + d*x)^3*(A + B*cos(c + d*x) + C*cos(c + d*x)^2))/(a + a*cos(c 
+ d*x))^(5/2),x)
 
output
int((cos(c + d*x)^3*(A + B*cos(c + d*x) + C*cos(c + d*x)^2))/(a + a*cos(c 
+ d*x))^(5/2), x)